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Abstract 

Starting from the QCD Lagrangian, we present the QCD Hamiltonian for 
near hght cone coordinates. We study the dynamics of the gluonic zero modes 
of this Hamiltonian. The strong couphng solutions serve as a basis for the 
complete problem. We discuss the importance of zero modes for the confine- 
ment mechanism. 
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I. INTRODUCTION 



The formulation of light front QCD is one of the most innovative enterprises in recent 
theoretical hadron physics. It resumes the pioneering efforts of the seventies in the parton 
model Its intention is to connect the successful parton model at large resolution 

with the constituent quark picture of hadrons appearing in spectroscopy. The new start |^ 
is not without a knowledge of the problems which have been experienced in the first works. 
It is well documented that renormalization in Hamiltonian field theories is, with currently 
available methods, more cumbersome than in covariant descriptions. A naive gauge fixing 
procedure on the light front leads to an easy resolution of Gauss' Law. However, this naive 
method is not correct - the correct method involves so-called 'zero mode' degrees of freedom 
dependent on the transverse coordinates. These zero modes cannot be gauged away, and 
become an integral part of the dynamics. In addition we expect that the nontrivial vacuum 
structure evident from equal time quantization brings new induced couplings into the light 
front Hamiltonian. 

Our investigation concentrates on the role of the zero mode fields. We start from the near 
light front frame advocated by the St. Petersburg and Erlangen groups and introduce 
a finite volume (L x L x L in the spatial directions), thereby controlling possible infrared 
singularities. This choice of coordinate system involves quantization on a space-like surface, 
which makes it easier to relate the occurring phenomena to equal time physics. However, 
all the usual complexities of negative energy states and nontrivial vacua are also present. 

An axial gauge is natural in the infinite momentum frame: the axis of motion singles out 
a preferred direction. The chosen gauge d-A_ = is actually a variant of the axial gauge, 
since the latter is incompatible with the boundary conditions. Consequences include the 
appearance and quantization of the zero-mode gauge fields a_ which depend only on the 
two transverse coordinates and can be chosen color diagonal. The asymmetrical dependence 
of the zero modes on only the transverse coordinates coincides with the asymmetry of the 
space coordinates in this near light front frame. Explicitly, the light front spatial x~ direction 
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stems from mixing the former spatial coordinate with the former time 

The zero modes degrees of freedom appear in a transverse Hamiltonian which is coupled 
to three-dimensional dynamics via the fermions, transverse gauge fields and the light front 
Coulomb law. In the strong coupling limit the kinetic term of the transverse Hamiltonian 
becomes dominant and, in this limit, the Hamiltonian is identical to a Hamiltonian describing 
independent rotators at each lattice site. This is the starting point of our investigation. Using 
basis functions according to this dominant kinetic term, we make an expansion of the zero 
mode Hamiltonian. In the strong coupling basis we also evaluate the Coulomb term and 
find indications that the interaction of external sources is confining. 

In Hamiltonian formulations quantized on the exact light front, the zero modes a^{x±) 
correspond to high energies. They are infinitely long wavelength excitations with respect 
to the spatial variable x~ , and, simultaneously, they are very short wavelength (ultraviolet) 
excitations with respect to the time variable x~^. Therefore some approaches [^] include 
the zero mode physics in their renormalization program, without explicitly solving their 
dynamics. In our approach of near light cone coordinates, the zero modes are independent 
degrees of freedom and retained after the solution of Gauss' Law, since they correspond to 
gauge invariant quantities (as evident in their role as eigenphases of the Polyakov loops). 
Resolution of the Gauss law constraint indeed does not permit the elimination of their 
conjugate momenta p~{x±). 

The existence of zero modes indicates that the local color charge of all external sources 
will be singlet - all hadronic bound states must be color singlets. This is the first require- 
ment of color confinement. The second requirement for color confinement dictates that the 
interaction energy increases with increasing separation. To make a reliable calculation in 
this case, one has to know, besides the light front Coulomb potential, how gluon fields prop- 
agate in the 'background' field a_{x±) over large distances. It is argued that on a strongly 
coupled lattice the background field modifies the behavior of transverse plane gluon waves 
with color charge. Due to the fluctuations of the background they are limited to propagate 
over short distances. Thus they cannot cancel the linear confinement potential induced by 
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the gauge choice. In this way, the axial gauge supports confinement from the beginning 
and therefore seems to be the best starting point for QCD. However, confinement is not so 
obvious for neutral transverse gluon fields, whose two-dimensional longitudinal zero mode 
part has been eliminated in the procedure of solving the Gauss law constraint. 

Light front gauge theories have been set up using the so-called near light front frame 
@J^J^,|], which enables one to study the approach to the exact light front. The tilted 
coordinates are defined as 



X = —pz [X — X 



The transverse components x^ and are unchanged; x~^ is the new time coordinate and x~ 
is the remaining spatial coordinate. As finite quantization volume we will take a torus and 
its extension in x~, as well as in x^, direction is L. The metric tensor reads 
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where /i , = +, 1, 2, — . It defines the scalar product of two 4-vectors x and y: 



:i.2) 



x^y^ = X y^ + x^y — rfx y — x±y± 



x_y+ + x+y_ + 7] x+y+ - x±y± 



;i.3) 



where 



x±y± = x^y^ + x^y"^ 



[lA) 



Obviously, the exact light-cone is approached as the parameter rf goes to zero. For non- 
zero 7], the transition to the near light front coordinates from an equal time frame can be 
formally identified as a Lorentz boost combined with a linear transformation which avoids 



time dependent boundary conditions ||^, as required in the canonical formulation. The boost 
parameter (3 for this Lorentz transformation (in the direction) is given by 

indicating that for 77 the relative velocity of the two frames v = (3c ^ c{= 1). This 
is connected to the interpretation of the near light front frame in terms of the infinite 
momentum frame. As discussed previously, the use of these near light front coordinates 
allows us to quantize the theory on a spacelike surface at equal light front time. This 
is important, as the ends of our 'box' are spacelike separated, and can thus exchange no 
information. Note the contrast to conventional light front coordinates, in which the ends of 
the 'box' are separated by a light-like interval which, in turn, implies that the surface of 
quantization contains points which can be causally connected. 

The spectrum of massless partons p^p^ = with transverse momentum p± can be easily 
calculated from their dispersion relation which reads in the near light front coordinates 

2p_p+ + 7]\p^y-pl = 0. (1.6) 

In our coordinates, p+ is the energy variable (being conjugate to x"*") and p_ is the longitu- 
dinal momentum variable. This dispersion relation has two p+ solutions, an upper branch 
p+,up and a lower branch 

P+,up/down = -^{P-T \/{P-y + V^Plj ■ (1-7) 

In the limit of small t], Eq. ( |1.7| ) goes over into the conventional light front energy 

P+,np P+,lc = pl/2p_,ic. (1.8) 

In Figure |I|, we show these two branches together with the light front Hamiltonian 
p+^ic = p]_/2p-^\c. We have taken the maximum transverse momentum p_L = 1/a, 7] = 0.1, 
and all energies and momenta are in units of 1/a, where a is the lattice spacing. Note that 
negative momentum states correspond uniquely to negative energy states in the conventional 



light front formalism, but with our choice of coordinates, this clean division no longer holds. 
However, it is true that for small r], the negative momentum states with positive energy are 
of very high energy, indeed. One possible procedure to find the momenta for regularization 
of these high energy states is to calculate the longitudinal momenta corresponding to a 
maximum absolute energy 



For p± = 1/a, this cutoff corresponds to p_ = in p+,up/down and renders the dispersion 
relation single-valued once again: a unique energy p+ corresponds to each momentum p_ . 

Since we are finally interested in using p+^ic as our effective Hamiltonian, the upper energy 
cutoff gives the minimal p_ momentum where this approximation still makes sense. When 
we use the expression in Eq. ( |1.8|) for p+,ic we obtain a minimal p_ = f]/2a, resulting from our 
energy cutoff. To introduce an efficient light front effective theory, it is necessary to eliminate 
partons with negative p_ states with special attention to the partons with < j9_ < ri/2a. 
This procedure will not be discussed here, but we believe that this elimination is very 
important in order to obtain a constituent quark picture on the light front. We would 
further comment that our effective Hamiltonian approach discussed in the next chapter may 
provide a workable framework to include these high energy modes. 

This paper is divided in six sections. In the next section we review the QCD Hamiltonian 
in the case of these tilted light front coordinates. In section 3 we give the solution to the zero 
mode sector in the strong and weak coupling approximations. A calculation of the ground 
state energy in the two site approximation is presented in section 4. Section 5 contains a 
discussion of confinement. Finally, we summarize in section 6. 

II. QCD HAMILTONIAN NEAR THE LIGHT FRONT 

Canonical formulations of QED in the axial gauge and in the light front gauge have been 
developed in analogous ways - starting from the respective canonical Weyl gauges. After the 
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implementation of the Gauss law constraints, the resulting Hamiltonians appear to be rather 
similar 0. Moreover, the QCD Hamiltonian in axial gauge representation has recently been 
derived [0. Here we will outline the derivation of the near light front QCD Hamiltonian, 



which has been given by the Erlangen group |TT| also. We restrict ourselves to the color 
gauge group SU{2) and dynamical gluons; only an external (fermionic) charge density pm 
is considered here. General SU{N) results, including dynamical fermions, are given in the 
Appendix. 

The Lagrangian in the near light front coordinate system reads 



— _)_ I F« P« _|_ !L F« P« I _ 1 pa pa _ a /.a 



(2.1) 



where the color index a is summed from 1 to 3, and the transverse coordinates are labeled 
hj i = 1,2. We will also use the matrix notation; for example A_ = A"r"/2, where the 
r" are the SU (2) matrices. The field strength tensor contains the commutator of the gauge 
fields: 

i^M- = d^A, - d,A^ - ig[A^, A,] , (2.2) 

with the coupling constant g. 

The A";^ coordinates have no momenta conjugate to them. As a consequence, the Weyl 
gauge A°^ = is the most natural starting point for a canonical formulation. The canonical 
momenta of the dynamical fields At-, A^ are given by 

djC 

- ~ dFl_ ~ +-' 
n- = ^ = F% + rt'Fl, (2.3) 

From this, we get the Weyl gauge Hamiltonian density 

Hw = ^n"_n"_ + If^2Fi2 + a E (n- - F^f . (2.4) 

We choose periodic boundary conditions in x~ and x± on intervals of size [0,L]. Using 
the appropriate periodic delta functions, the quantization is straightforward. However, the 



Hamiltonian has to be supplemented by the original Euler-Lagrange equation for as 
constraints on the physical states 



= (D^^n^ 1$) = 0, (2.5) 

with the covariant derivatives 

D"^^ = + gr'''A1, 

Df = d^6^^ + gT'^Al, (2.6) 

where /"'^'' are the structure constants of SU{2). 

These equations are known as Gauss' Law constraints. Since the Gauss' Law operator 
commutes with the Hamiltonian 

[G''{x^,x-),Hw]=0, (2.7) 

time evolution leaves the system in the space of physical states. Furthermore, H^r is invariant 
under time independent residual gauge transformations whose generator is closely connected 
to Gauss' Law [^]. 

In order to obtain a Hamiltonian formulated in terms of unconstrained variables, thus 
rendered available for approximations without breaking local gauge invariance, one needs 
to resolve the Gauss' Law constraint. This is an important step, since it frees us from 
implementing additional restrictions on operators or states. A Hamiltonian expressed in 
terms of unconstrained variables appears more complicated, but it is our belief that there is 
much to gain from this apparent increase in complexity. 

Via unitary gauge fixing transformations [p|,p!0| one indeed can achieve this resolution 
with respect to components of the chromo-electric field. These transformations render a 
Hamiltonian independent of the conjugate gauge fields. In other words, the latter become 
cyclic variables. Let us choose the ' — ' (minus) components as the variables to be eliminated. 
Classically this would correspond to the light front gauge A_ = 0. However, this choice is 
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not legitimate for our setup in a finite box. Only the (classical) Coulomb light front gauge 
{d-A_ = 0) is compatible with gauge invariance and periodic boundary conditions. The 
reason is that carries information on gauge invariant quantities, such as the eigenvalues 
of the spatial Polyakov (Wilson) loop 



V{x±) = Pexp 



(2.8) 



ig J dx A_(x_L,x ) 

which can be written in terms of a diagonal matrix a-{x±) 

p(x±) = 1/exp [igLa_{x±)] . (2.9) 

Thus, we obviously need to keep these 'zero modes' a-{x±) as dynamical variables, while the 
other components of A_ are eliminated. The zero mode degrees of freedom are independent 
of x~ and, therefore, correspond to quantities with zero longitudinal momentum. 

The unitary gauge fixing transformation can indeed be chosen in such a way that A_ 



becomes cyclic, apart from the zero modes mentioned |TT|. In order to eliminate the conju- 
gate momentum, n_, by means of Gauss' Law, one needs to 'invert' the covariant derivative 
D_. After the unitary transformation D_ simplifies significantly (compare to Eq. ( |2.(j| )) 

^ d- = d- -ig[a-, , (2.10) 

where this is understood as an operator equation: e.g. -D_/ —>■ d^f = d-f — ig [a_, /] for 
arbitrary f. Now Gauss' Law can be readily resolved: in the space of physical states one can 
make the replacement 

U_ix^,x.) ^ p.ix^) - [dz') G^{x^,y-). (2.11) 



The inversion of d^ can be explicitly constructed in terms of its eigenfunctions, cf. [[TO| , p 
The zero mode operator P-{x±) is also diagonal and p^{x±) is the conjugate momentum 
to a^{x±). It has eigenvalue zero with respect to i.e. d-P- = 0, and is therefore not 
constrained. 

The appearance of the zero modes also implies residual Gauss' Law constraints. In the 
space of transformed physical states they read 



/ dx-Gl \x) = / dx- {Dfu\_ + gpl) \x) = 0. (2.12) 

These two-dimensional constraints can be handled in full analogy to QED, since they corre- 
spond to the diagonal part of color space. This further gauge fixing in the SU (2) 3-direction 
is done via another gauge fixing transformation, which leads to the Coulomb gauge repre- 
sentation in the transverse plane for the neutral fields. In other words, we eliminate the 
color neutral, a; "-independent, two-dimensional longitudinal gauge fields 

1 f 

= j^j dy-dy^d{x^ - yj_)Vi_ (V^ • Ai{y^, y_)) - . (2.13) 

Here we use the periodic Greens function of the two dimensional Laplace operator 

d{z^) = E ^e^'"'" ' Pn = , (2.14) 

where n = (ni,n2) and ni,n2 are integers. The conjugate momenta of these fields, p^^{x±), 
are defined analogously. Resolution of the residual Gauss' Law allows one to replace them, 
in the sector of the transformed physical space |$'), by the neutral chromo-electric field 

e^{x^) = gV^ J dy-dy^d{x^-y^) [f'^'Al{y^,y-)Ii'^{y^,y-)+pl{y^,y-)Y^. (2.15) 

At this point it is convenient to introduce the unconstrained gauge fields and their conjugate 
momenta: 

A'^{x^, x~) = A^{x^, x~) - ai(x_L), 

n'^(xx, X-) = nx(xx, X-) - -^pi(xx). (2.16) 

These relations turn out to be important for neutral gluon exchange; recall that the sub- 
tracted fields are diagonal in color space. Note that the physical degrees of freedom and 
n'j^ still contain (x_l, a;~)-independent, color neutral, modes. Therefore, there is a remnant 
of the local Gauss' Law constraints - the global condition 

Q^W)^ I dy-dy^{f'^'Al(y^,y-)U',(y^,y-)+pl(y^,y-)}\^')^0. (2.17) 
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The physical meaning of this equation is that the neutral component of the total color 
charge, including external matter as well as gluonic contributions, must vanish in the sector 
of physical states. 

The final Hamiltonian density in the physical sector explicitly reads 

n = tT [^^A', - d^A, - ig[^,, A^]f + Itr [H'^ - {d^A!^ - ig[a.,A^])f 



1 

+ ^tr 



2 



+ -;^P-\xA_)pi{xA_) 



where p and q are matrix labels for rows and columns, a_g = {aS)qq and the prime indicates 
that the summation is restricted to n 7^ if p = g. The operator G'^ {x±,x~) is defined as 

G'^ = V^n'^ + gr'^'^A' ^ (n' 1 - ie l) + gp.^ . (2.19) 

In Table 1 we show the number of degrees of freedom for the different stages of the 
gauge fixing procedure. In order to count, we discretize configuration space as sites. 
Constraints of course reduce the number of independent variables and are therefore sub- 
tracted. Originally there are three vector components of the gluon field in three colors at 
each site. Gauss' law represents three (color) constraints at each site. After the first unitary 
transformation we arrive at a intermediate Hamiltonian which is not explicitly given in the 
text. The number of degrees of freedom as well as the number of constraints have been 
reduced. Note that the residual Gauss law is two-dimensional and color neutral. Finally, 
we arrive at a formulation where there is only one global constraint left. Concomitantly, the 
final A'j_ still contains the global color neutral zero mode which has not been subtracted (cf. 
Eqs. ( p.l3| , |2.14D ). Thus, one explicitly sees that the number of independent unconstrained 



degrees of freedom is 6A^'^ at any stage of the formal development. 
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TABLES 

TABLE I. Counting the degrees of freedom in the SU{2) Hamiltonians 



Original Ti. 






3 * 3 * Ar3 


Gauss' law constraints: 


-3*N^ 


Intermediate H, 




74f, Af(x_L,x_) : 




Zero modes a^_{xx) : 




Residual Gauss' law constraints: 




Final H 




A'f,A'2^{x^,x_) : 


2*3* - {N^ -I) 


Zero modes a^{x±) : 




Global neutrality constraint: 


-1 
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The formulation of axial gauge near light front QCD for SU{2) gauge fields is complete 
at this point. For details of the rather subtle derivation we refer to original references 



|Ty,[Tl]|. See also ||T2| for a pedagogical review and ||13| for a critical discussion focusing on 
topological aspects. In [|T^ the method using unitary transformations has been applied to 
QED quantized on the exact light cone. The above Hamiltonian already is rather complex 
and nonlocal, just as the most familiar example of a gauge fixed theory, Coulomb gauge 
QED. Despite the complexity it serves as a promising starting point for further studies since 
approximations can be made without breaking local gauge invariance. 

III. ZERO MODE DYNAMICS 

The principal advantage of an exact light front formulation is the apparent triviality 
of the ground state which simplifies calculations of the hadron spectrum. The light front 
vacuum, however, is not guaranteed to be trivial in the zero mode sector. In using the near 
light front coordinate system, we can study the complex zero mode structure influencing the 
dynamics of long distances. The zero mode sector in light front physics differs from equal 
time Hamiltonian physics where the long range physics is low energy physics. As can be 
seen from the dispersion relation for massless particles on the light front, = soft 
modes (states with small momenta) become high energy states. In this way, high energy 
physics becomes tied to long range physics, contrary to the equal time formulation. This 
physics appears in deep inelastic scattering at small scaling variable and is related to the 
long distance features of the proton. We will focus on the zero mode sector in order to try 
to acquire some insight into its dynamics. Therewith we hope to obtain a sound basis for 
further numerical studies. 

From the comparison of abelian and non-abelian theories, striking differences show up 
in the zero mode sector. Recently, in the equal time formalism, the zero mode sector in 
QCD has been claimed to be relevant for the confinement phenomenon |T3. On the level of 



approximations and restrictions followed below, the formal differences between light front 
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and equal time approach are rather small and, consequently, results and methods are similar. 

In this work we do not restrict ourselves to the strong coupling approximation. We will, 
however, start with the strongly coupled theory to define our set of basis functions. As 
before, we will restrict ourselves to SU (2) without matter - pure gluonic Yang-Mills theory. 
It already has the typical non-abelian features such as the Coulomb term which explicitly 
contains the zero modes in the denominator and the non-standard kinetic energy for the 
zero modes. 

The zero mode degrees of freedom couple to the three-dimensional gluon fields via the 
second term in Ti. shifting the longitudinal momenta of the transverse gluon fields (Eq. 
( p.l8| )). They affect the Coulomb term and the two-dimensional electric fields e±. The 
latter coupling is typical for the light front and is absent in the equal time case. We neglect 
these couplings and consider the pure zero mode Hamiltonian 



h = d X 



(3.1) 



This Hamiltonian is obtained from Eq. ( p.l8| ) with only the zero mode and its conjugate 
momentum retained, and an integration over the longitudinal variable. We recognize 'elec- 
tric' and 'magnetic' contributions in h, the zero mode Hamiltonian - the first and second 
term, respectively. The light front variables mix the ordinary spatial and time variables so 
the labeling above is to be understood in analogy with the equal time Hamiltonian. 

Even at this level of approximation, this zero mode Hamiltonian differs from the cor- 
responding one in QED. The reason is the hermiticity defect of the canonical momentum: 
pL 7^ • This might seem to be a strange property for a momentum operator, but is per- 
fectly allowable, in analogy with the naive radial Schrodinger momentum operator which is 
also non-Hermit ian. 

For notational simplicity, we now omit the color index and work with the Schrodinger 
representation of Eq. (|3.1|) 



h 



J d'^x± 



(3.2) 
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where J {a J) is the Jacobian and equals the Haar measure of SU{2) 



J{a^{x±)) = sin^(^a_(x_L)). 



(3.3) 



The Jacobian is connected to the hermiticity defect of p_ ; they stem from the gauge fixing 
procedure taking into account the curvihnear coordinates. The measure also appears in 
the integration volume element for calculating matrix elements. For ease of calculation, we 
introduce dimensionless variables 



in which Hamiltonian and Jacobian respectively read 
g^L 1 S ^, , 5 



h 



J (fx± 



8 Jiip{x±)) 5ip{x±) 



+ 



6(p{x±) rfg'^L 



(3.4) 



(3.5) 



(3.6) 



As in earlier approaches, see e.g. [|T7|, (p will be treated as a compact variable, < ip < tv. 
At this stage it is necessary to appeal to the physics of the infinite momentum frame 



to factorize the reduced true energy h^ed and the Lorentz boost factor ^ 



2r?' 



smce 



essentially h is a light front energy, and it is well known how these behave under a Lorentz 
transformation. We rewrite 



with 



h 



red 



J cfx± 



2ri 



(3.7) 



Ji- + 



-(Vx^)^ 



(3. 



4 J Sip 5ip g^Lrj 

It should be noted that the coefficients of the two terms are reciprocals of each other. 

Since the integral over transverse coordinates can contain arbitrarily small wavelengths, 
we have to regularize the above Hamiltonian hj-ed- We do this by introducing a lattice 
to evaluate the transverse integral. The lattice vector b numbers the lattice sites, and ei 
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and 82 are the two unit vectors on the two-dimensional lattice. In order to have standard 
commutation relations on the lattice the derivative on the lattice becomes = t-t — ro^. 
We further explicitly pull out the dependence on the lattice cutoff by defining a new reduced 
Hamiltonian /ij-cd and substituting t] = -^aM, where M is a typical hadronic mass (see, e.g. 

my- 

h = ^hred , (3.9) 
Irja 



with 



Since the effective coupling constant, 

9e. = ^, (3.11) 

contains the large factor LM, the product of lattice size in the longitudinal direction and 
the hadron mass, a strong coupling approach seems to be a good starting point. 
We separate the Hamiltonian into electric and magnetic contributions: 

hrea = T.H + T.hf, (3.12) 

b b 

where 

/^g=-^eff7T^^T^, (3.13) 

and 

hf = ^UVb-^b^f- (3.14) 



Note that we do not introduce 'radial wave functions' nor effective potentials as in |T0|,|T5 
For each lattice site 6, the electric part of the Hamiltonian /i2 (the kinetic energy) has the 
Gegenbauer polynomials C„g (</?/;) as eigenfunctions: 

hlCn,{^i) = glsn^n^, + 2)Cn,{n), (3.15) 
16 



with 

C„,(,,).^|!MK±iM|, (3,6) 

and 

J{(p)Cni(p)Cm{(p)d(p = Snm ■ (3.17) 

The strong couphng wave functions of the full transverse lattice are product states charac- 
terized by a set of quantum numbers {n} = {rig}, 

^{n}i^)=l[Cn,i^b)- (3-18) 
b 

These functions form a complete and orthonormal basis for the zero mode sector. They 
satisfy the energy eigenvalue equation 

J2hl^{„,yi^) = g'^sE^b(^b + '^)'^{n}iv)- (3-19) 

b b 

The ground state in this limit corresponds to all = - a constant wave function 

^{0}(^) = n\/^' (3-20) 
b * 

and the ground state energy is zero 

Eo = 0. (3.21) 
The first excited energy level is A^^-fold degenerate - an excitation at a single lattice point 



In strong coupling this level is separated by a large amount from the ground state energy 

El = ?>gl^. (3.23) 



So far our results are equivalent to those of []T3| to within re-definitions of wave functions and 



integration measures. In |jT^ weak coupling variational solutions for the full SU{2) lattice 
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Hamiltonian are given. Furthermore, studies in (l+l)-dimensional Yang-Mills theory |T8 
give formal extensions to construct wave functions for SU{N) gauge theories. 

The magnetic term of the Hamiltonian couples nearest neighbor lattice points. In the 
strong coupling limit its contribution may be obtained perturbatively (as it has the coefficient 
l/^fgg) by evaluating it with the basis function of the ground state. The result of this is 

(^{0}l E hfl^w) = 7^ - iV • (3.24) 
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Since this energy is proportional to A^^ = {L/af' and g^^^ grows linearly with L, this part 
of the zero mode dynamics represents a negligible surface effect for the three-dimensional 
system in the strong coupling approximation. 

Next, we discuss the weak coupling limit — 0. In this case we can simplify the kinetic 
term of the Hamiltonian by defining new variables ap 

= — , (3.25) 
Kg 

with = \/2LM . Then the reduced Hamiltonian becomes: 

ke. = E -?(^ --'(-^«.^)^ + E(«^(^^) - -^Uf- (3-26) 
Expanding this for small g, we obtain 

+ (3.27) 

The eigensolutions of this Hamiltonian are known to be spin waves. Going over to Fourier 
momentum representation, 

ag = Ee^''%, (3.28) 

k 

with ki = 27Tni/Nj_a and = 0, ±1, ±2, ... we have 

- ^ {- (|| + + 4E.i.>^ f . (3.29) 

The eigensolutions ipK of /ired in the weak coupling approximation are decoupled harmonic 
oscillators for each k, with frequencies 



Ijj 



he 



4Esin'y. (3.30) 



Because of the 'radial Laplacian' it looks as if the eigenf unctions would have to vanish at 
the origin to be normalizable. However, as in the Schrodinger equation in three dimensions, 
the Jacobian J allows a constant wave function at the origin. Consequently, the eigenvalue 
of is given by the sum over the modes: 



k 



\ 



4i:sm^|f|. (3.31) 



which gives in the A^_l oo limit spin waves with cug = y /ci + fc^ . 

In the weak coupling limit the zero mode Hamiltonian supports solutions similar to 
QED. The strong coupling limit, however, yields different results: 'gluonic' excitations are 
suppressed because of large energy gaps. This is due to the Jacobian, which can be traced 
back to non-abelian self interactions in the original Lagrangian. 



IV. TWO SITE TRUNCATION 

We have obtained solutions in both the weak and strong coupling regime. We will now 
study the problem in the intermediate region. As an initial effort, we will not solve the 
problem for the full lattice. Rather, we will calculate with what is essentially a cluster 
expansion [|T^, and we will start with the simplest, two-site cluster, in which either site 
(or both) can be excited to high energy states. We will obtain the solution for the low- 
lying spectra of the system approximated as a low density of excitable two-site clusters. 
This method can be envisaged as the starting point of a more ambitious Hamiltonian based 
renormalization group technique, like the contractor renormalization group method (CORE) 
li, or |0 . 



We handle the calculation of the energies via an effective Hamiltonian method. We work 
in the representation of the strong coupling solution of /i| and divide the two-site subspace 
into a P and Q space, such that P + Q = 1 with 
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^ = {|o,o)}, 

Q = {\n,m); n,m^ 0,0}, (4.1) 

where n, m represent the indices of the Gegenbauer polynomials. Note that we have picked 
our P space as the strong coupling two-site ground state. Then the two-site energy E2 is 
given by the non-perturbative solution of the Hamiltonian in Eq. ( p.l2|) , truncated to two 
lattice sites. Explicitly, this Hamiltonian is: 

h2 = h' + (4.2) 



with 



and 



/ /1J_J J_ + U_j J_| (4 3) 



= ^(^1 - ^2?, (4.4) 

9eS 



where the subscripts label the sites. 

Within the effective Hamiltonian method, the two-site energy is given by (see, for ex- 
ample [PH] ): 



E2 = Ph2P + Ph2Q — Qh2P. (4.5) 

E2 - Qh2Q 

The self-consistent solutions of this equation provide the low-lying spectra in this method. 
The strong coupling basis states are eigenstates of h'^: 

h^\n, m) = {n{n + 2) + m(m + 2)} \n, m) . (4-6) 

Thus, the non-trivial matrix elements are those of /i™, and are of the form 

{n,m\{ipi- Lp2f\n ,m!) . (4.7) 

Although our states are two-site states, the operators appearing in the matrix elements 
are simple one-site operators, and thus we can consider the states to be products of one-site 
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states. This reduces the evaluation to sums and products of one-site matrix elements, which 
are given as: 

{n\'^\n') = — ioT n = n' 



2 1 1 



{n\ip\n) 



for n + n' = odd 



TT V("+"'+2)2 (n-n')2^ 

for n + n' = even, n ^ n' 



2 



/ I 2| /\ I TT I D [2(n+l)l-j / A o\ 

{n\Lp \n) = \ " (4.8) 

Before performing any numerical calculation with the effective Hamiltonian, it will be illu- 
minating to investigate the strong and weak coupling limits of this theory. 

In the strong coupling limit we expect to obtain the result of perturbation theory in 
l/gcs- The energy in this limit is easily calculated: 

E, = (0, 0|/z-|0, 0) = ^ - 2 {{Om?} = ^ f^^ - l) • (4-9) 

9eS ycff V " / 

In the weak coupling limit we can solve the Schrodinger equation for two neighboring sites 
- the Hamiltonian will simply be the two-site version of the earlier spin wave Hamiltonian, 
Eq. ( |3.27| ). We call the respective variables ag^ = x and a^^^ = y , then we have to find the 



for n = n' 



eigenenergies of the Hamiltonian: 

h =-( — + -(— + -— \+{x-yf (4 10) 

'^'^ \ dx"^ X dx J \ dy"^ y dy J 

It should be noted that this Hamiltonian is invariant under x ^ y, and thus the eigenfunc- 

tions \l/2(x,y) can be chosen to be symmetric under the interchange of x and y (\l/2s(x, y)), 

or antisymmetric under the interchange of x and y {'^2a{x,y))- Each of these symmetric or 

antisymmetric sets of solutions form a 'tower' of excitations. The first symmetric excited 

state becomes degenerate with the ground state of the original problem in the weak coupling 

limit, and the first antisymmetric state has a greater energy than the first symmetric state. 

As usual one factorizes the wave function 

^2{x,y) = —<l>2{x,y) , (4.11) 
xy 
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resulting in the Schrodinger equation 

The center-of-mass motion is then separated: 

$2(x,y)=e^^\2(r), (4.13) 

with R = {x + y)/2 and r = x — y. The Hamiltonian corresponding to the relative motion 
(r) is a simple radial harmonic oscillator: 

W. = -2|^ + rl (4.14) 

The lowest states of the symmetric and antisymmetric 'towers' are solutions to this Hamil- 
tonian. The energies of these states can be read directly from Eq. ( |4.14| ); = V^, and 
E2a = 3-\/2, respectively, giving a energy gap between the states of 2\/2. 

Thus, the results for the energy gaps of the low-lying states in the weak coupling limit 

are 

E2s ~ -^'ground = , 

^2a- Aground = 2^2. (4.15) 

We now proceed to calculate the low-lying spectra via the effective Hamiltonian method. 
In the numerical calculations, we cannot keep all states in the Q space - our choice is to 
cut at a high two-site energy, calculate E2, then increase the size of the Q space to check 
for convergent results. This procedure was carried out for each choice of coupling constant 
ges, and the typical number of two-site states kept in the Q space at convergence was about 
300. 

The numerical solution of Eq. ( |4.5| ) for E2 is given in Figure |^. In the strong coupling 



limit (large g'^g) the large gaps in energy are evident, and the numbers agree with the 
unperturbed energy of the states, given in Eq. ( [4.61 ). In this same limit, the slope of the 
ground state energy as a function of the inverse square coupling agrees with the analytic 
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calculation of Eq. ( [4.9|) . In the weak coupling limit, the results for the gap energies were 
Richardson extrapolated for the l/(7cfr ^ ^ limit. This extrapolation matched the analytical 
results of Eq. ( [4.15[ ) to five significant figures. Thus, we have obtained two-site solutions 
for the entire range of coupling which agreed with analytic results in the weak and strong 
coupling limits. 

It is remarkable that we succeeded a one-dimensional strong coupling basis throughout 
the range of coupling strengths. Results for the spectra of A^-sites are straightforwardly 
obtained as long as the number of excited two-site clusters is small compared with N/2. 
This is the 'low density' approximation. 



V. EFFECT OF ZERO MODES ON CONFINEMENT 

In order to study the confinement problem, we start with the assumption that the three- 
dimensional gauge Hamiltonian can be treated in a weak coupling approximation with small 
gauge coupling gf^, whereas the physics of the two-dimensional zero mode subsystem can 
be obtained in a strong coupling approximation for g^^ = ^ ■ We will see that the 
asymmetric treatment of transverse and longitudinal spatial coordinates has effects, which 
go beyond the violation of rotational symmetry in strong coupling lattice gauge theory. They 
are connected with the procedure of initially choosing an axial gauge in three dimensions 
and then letting a two-dimensional Coulomb gauge follow. One advantage of this approach 
is that the dynamical role played by the zero modes becomes particularly illuminating. The 
zero modes preserve the axial linear confinement in first order perturbation theory. Thus, 
we argue that already at the level of the zero mode sector of SU (2) pure glue QCD we see 
evidence for confinement. 

Let us demonstrate the usefulness of the basis functions obtained in the strong coupling 
approximation, cf. Eq. ( |3.16| ), by considering matrix elements of the Coulomb term. In 



terms of the variables defined on a discretized transverse space the Coulomb potential 
explicitly reads 
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27rn(z —y )/L 



7m + {p- q)^i 



(5.1) 



In order to separate possible singularities we work out the p, q sum: 



2 -y 



mm 



4 —JO 

b 







+ ^T. Tdz- r dy J2 

4 "T^JO JO „ 



'^±21 Z')G'±i2ib,y') i2TTn{z-^y-)/L 



(5 



vrn - 

Obviously the 'abelian' terms with G[i and (j22 ^ire infrared regular since tt, 7^ in that case. 
The non-abelian terms can have singularities for ip —>■ 0,tc. Exploiting the strong coupling 
basis, we show that these terms are also infrared regular. Such a 'dynamical regularization' 
is anticipated, because of the connection of the non-standard zero mode kinetic energy, the 
hermiticity defect, and the Coulomb term: the Jacobian vanishes at the points, where the 
'propagator' becomes singular. Consider the last two terms in Eq. ( |5.2| ), which can be added 
using z ^ y, n ^ —n. This yields the non-abelian Coulomb contribution He. 



G1ji(S,j-)G'j 



l.l2\"^y ) ^2mn(z--y')/L 



vrn 



(5.3) 



, - - n V - ^bY 

In this expression, we can identify the 'Coulomb propagator' D^^z" — y^,ip^) in position 



space, and it can be evaluated with the result 



Dciz -y , <^5) = - E 



^2iTin(z ~y )/L 



2 V (™ - VbY 



^2i(pi-{z -y )/L 



2 sin" 



i{z -y )cotv3g 



(5.4) 



In the continuum limit one finds additional terms besides the linear propagator in one 
dimension. In the strong coupling approximation we integrate this Coulomb propagator and 
the off-diagonal Gj_2i^j_i2 with the Gegenbauer polynomials Co(v2) over dipJ{ip) appropriate 
for the curvilinear coordinates. At this point it is evident that the Jacobian indeed prevents 
possible infrared singularities mentioned above. 



24 



Let us discuss the first term Di in the Coulomb propagator. With the strong couphng 
ground state ^{0} it leads to a Coulomb energy 

rL rL 



(vl/|o}|i/«|^{o}) = (^{o}|a2E I I dy- G'^,,{b,z-)G'^,,{b,y-)D,\^{,}) 

b 

= La'Y: tdz- ['^dy~G'^,,{b,z-)G'^,,{b,y-) x 
^ JO Jo 



b 

T 

g27rj(2 ~y )/L _ ^ 



for large L, with 



L rL 

dz~ I 
Jo 

L 

2TTi{z~ — y 

La'Y.\Qi2m\ (5.5) 



Qi2{h) = rdz~G'^,,{b,z-). (5.6) 
Jo 



One sees that the expectation value of the off-diagonal 'charge' Quib) at each transverse 
site b has to vanish in order to avoid an infinite Coulomb energy in the continuum limit 

{Q,2(b)) = V6. (5.7) 

Recall that, in order to fully resolve Gauss' Law, the third component of the global charge 
should vanish in the physical sector (cf. Eq. (|2.171 )). These two conditions together suggest 



that the physical states are real color singlets not merely color neutral states with a color 
three projection equal zero. The second term D2 of the Coulomb propagator gives a linearly 
rising potential in the L ^ 00 limit. 



b 

E f^^dz- f^^dy- G'^2,{%z-)G'^,,{b,y-)\z- -y-\. (5.8) 
^ Jo Jo 



Neglecting terms proportional to g with gluons in G'^21 only considering the external 
charges one gets a confining linear potential. Two color spin 1/2 point charges coupled to a 
color singlet at the same transverse site interact in longitudinal direction with the potential: 

" ^^(^i^2)k" = ^g^\z~ -y-\^. (5.9) 
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The scale of the string tension is given in strong couphng by the lattice size a of the transverse 
lattice. In the continuum calculation it should be replaced by a correlation length generated 
in the transverse zero mode dynamics. The third term in the Coulomb propagator D3 does 
not contribute to the strong coupling ground state energy in the L ^ 00 limit, since then 

(gi2(&)) = 0. 

One of the main advantages of the light front Coulomb gauge or axial gauges for QCD 
is visible here. Whereas in QED the Coulomb gauge is designed to give the 1/r potential, 
the above gauge choices give confining potentials in zeroth order. Naturally, these potentials 
are linked to the choice of gauge, so one has to consider the perturbative corrections to the 
gauge potential. In axial gauge QED the potential for opposite charges has the form 

Mr^ = /|.--y-|l, (5.10) 

which appears to have the same confining properties as the QCD result. However, the first- 
order one photon exchange correction for soft photons with momenta {q-,qx) has a spin 
independent contribution 

i^^QED ^2 gj ^(g-) ~1 

M2 =g —, — . (5.11) 

Q-q- q- qi/q- 

After transforming M^^^ to coordinate space, it is seen that this one-photon exchange 
cancels the confining gauge artifact and the Coulomb potential 1/r plus spin dependent 
corrections remain. 

In QCD the onc-ghion exchange contributions have to be discussed separately for color 
charged gluons A'l^'^{x~ ,x±) and color neutral gluons. The A'}^'^ components of the transverse 
gluon fields interact with the zero mode a^{x±) fields. These 'background fields' have non- 

/I 2 

zero expectation values in the strong coupling limit for g^s- In this way the Aj_' fields may 
acquire a mass. In other words, we argue that the dispersion relation of transverse colored 
gluons possibly is changed to 



(-k. + ^k'_+ rj^ikl + ^{^')/L^)] . (5.12) 
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Consequently, the one-gluon exchange can no longer cancel the linear confinement potential 
at large distances, as the fields are now massive and thus finite-range. 

Our approach produces similar results as the similarity transformation scheme P,p3 



This scheme has been proposed by Glazek and Wilson for the light front Hamiltonian 
and by Wegner independently for condensed matter physics. It avoids vanishing energy 
denominators for the g_ — > region for one-gluon exchange by a cutoff A. In higher orders 
the region of validity for g_ can be enlarged successively (the cutoff A can be made smaller). 
In our approach near the light front we can give the physical origin of the cutoff A. It 
lies in the zero mode fields a-{x±) which limit the long-range propagation of gluon fields. 
We only use the strong coupling approximation for the ground state wave functional of 
the transverse zero mode lattice. Certainly higher orders and a more accurate description 
of the zero mode dynamics are necessary to prove that the linear confinement potential is 
preserved. Color neutral gluon fields A^(x~ , x±) do not directly interact with a^(x±), but 
by their construction in Eq. ( |2.16|) their two-dimensional longitudinal zero modes (g_ = 0) 
have been subtracted. The corresponding two-dimensional transverse parts, however, are 
still present as dynamical modes. 

It is more complicated to analyze two color charges separated in transverse direction. 
Consider two color charges separated by x± and oriented along the color 3-direction. As 
color 3-charges they experience the two-dimensional Coulomb potential via the electric field 

elix±) = gV± J dy~dy±d{x± - y±)pl,{y' ,y±). (5.13) 

The above propagator generates a logarithmic potential, which is also confining. This poten- 
tial is unaffected by the zero modes a_, since the strong coupling ground state expectation 
value of V_La_ vanishes. Explicitly, we obtain 

It corresponds to spreading flux lines in 2 dimensions and strongly violates the rotational 
invariance in the spatial coordinates. To investigate this problem further, one must under- 
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stand the role played by the zero modes in the continuum limit. Work in this direction is in 
progress |^ . 



VI. SUMMARY 

We have presented a near light front description of QCD. In our Hamiltonian approach, 
formulated in a finite volume, the modified light-cone gauge d-A^ = is a natural choice. 
The resulting asymmetric treatment of transverse and longitudinal dynamics matches the 
physics of deep inelastic scattering where large momenta of the hadron are involved. We 
argue that the well-known physical appeal of an infinite momentum frame can indeed be 
realized rather naturally in near light front coordinates. 

In the chosen canonical formulation. Gauss' law appears as a quantum mechanical con- 
straint to be implemented on physical states. It reflects the presence of redundant gauge 
variables. Using unitary gauge fixing transformations, these are eliminated by implementing 
Gauss' law. The result is a Hamiltonian formulated in terms of unconstrained dynamical 
variables. 

This Hamiltonian contains zero mode degrees of freedom, which only depend on trans- 
verse coordinates. The zero-mode dynamics apparently generates the dominant non- 
perturbative physics with interesting implications. In the strong coupling limit, g'^g = 
^ ^ 1, the energy gaps between excited 'gluonic' states and the ground state are large - 
in contrast to QED. Another indication of confinement is a linear quark-quark potential in 
longitudinal direction, which is not merely a gauge artifact - as it is in QED. In transverse 
direction confinement is supported via a logarithmic quark-quark potential. One has yet to 
discover how the two-dimensional dynamics relates to the full-three dimensional dynamics. 

The eigenfunctions of the zero mode Hamiltonian in the weak coupling limit are spin 
waves. We also have investigated the intermediate coupling regime by means of a cluster 
expansion. Via an effective Hamiltonian method the two-site energy was explicitly calculated 
for the entire range of couplings. We demonstrated that the results indeed converge to the 
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analytical expressions for g'^^ ^ 1 and g'^g <^ 1 . The spectra of A^-sites can be constructed 
in the 'low density' approximation. 

Our approach opens a potential path to construct constituent quark models on the light 
front by eliminating the negative energy solutions from the near light front Hamiltonian. The 
effective Hamiltonian method utilized here for the zero modes may also be the appropriate 
tool to study the mechanism of chiral symmetry breaking. 



APPENDIX: SU{N) HAMILTONIAN INCLUDING FERMIONS 

The resolution of the Gauss law constraint via 'unitary gauge fixing transformations' 



can also be achieved in the SU{N)-case, including dynamical fermions |p!0| , [Tl| . The final 
Hamiltonian in near light front coordinates is formulated in terms of the following SU (N) 
variables: 

• Fermion fields ip and ip^ obeying standard anti-commutation relations. (The N color 
indices, as well as Dirac and flavor labels are suppressed.) It should be emphasized 
that ip^ = |(1 — <y^)4' is a dynamical fleld, whereas on the exact light-front it is 
constrained. 

• Transverse gluon flelds, A'^^ii = 1,2; a = 1,....,A^^ — 1), without the neutral, two- 
dimensional, longitudinal parts, which have been eliminated. This means that 

di dx-A'i^'ixx^x-) =0, (A.l) 
Jo 

where the color index with subindex '0' refers to the Cartan subalgebra. In the usual 
representation this corresponds to diagonal A-matrices. The conjugate chromoelectric 
flelds are denoted with H^". 

• Zero mode gluon flelds a_ and their canonical momenta p_. Neither of these flelds 
depend on and both are color diagonal. Thus one can write 
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N-1 
co=l 

N-l ^co 

P-{x±) = J2 P-i^±)^ ■ 

co=l ^ 



(A.2) 



The Hamiltonian density in the physical Hilbert space exphcitly reads 

2 1 1 

H — (9_ — iga^) ■0_ — i-tp'^a^ ( V_l — igA'A i/j + m-ip^l3ip 



+ tr [d,A', - d,A[ - ig[A[, A',]f + -tr [nl - (d^A'^ - ig[a^,A]) 



1 

+ ^tr 



1 2 



-e_L — V_La_ 



+ 



2L2 



Co 



1 /-i _ 



^i27rn(z — j/ )/L 



(A.3) 



with p,q — 1, ....,N. The operator G'j^, the neutral chromo-electric fields e±, and the neutral 
charge Q'^", which appears below, are defined as the generalizations of the corresponding 
quantities in the main text. Note that the matter density contains the dynamical fermion 
fields, 



A" 



(A.4) 



The variables defined above are not constrained apart from the global conditions 



(A.5) 



The neutral components of the total color-charge vanish in the sector of (transformed) 
physical states. This completes the SU{N) formulation of axial gauge light cone QCD. 
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FIG. 1. Dispersion relations in the near light front coordinate system and directly on the light 
front. 
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FIG. 2. Energies of first 3 states in the zero mode sector calculated via effective Hamiltonian 
method. 
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